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Abstract 

This paper is devoted to study the gravitational charged perfect 
fluid collapse in the Friedmann universe models with cosmological con- 
stant. For this purpose, we assume that the electromagnetic field is so 
weak that it does not introduce any distortion into the geometry of the 
spacetime. The results obtained from the junction conditions between 
the Friedmann and the Reissner-Nordstrom de-Sitter spacetimes are 
used to solve the field equations. Further, the singularity structure 
and mass effects of the collapsing system on time difference between 
the formation of apparent horizons and singularity have been studied. 
This analysis provides the validity of Cosmic Censorship Hypothesis. 
It is found that the electric field affects the area of apparent horizons 
and their time of formation. 

Keywords: Electric field; Gravitational collapse; Cosmological constant; 
Friedmann models. 
PACS: 04.20.-q; 04.40.Dg; 97.10. CV 



*msharif@math.pu.edu.pk 
t abbasg9 1 @yahoo . com 



1 



1 Introduction 



Gravitational collapse of a massive star is the result of its self gravity. It oc- 
curs when the internal nuclear fuel of the star fails to supply sufficiently high 
pressure to counter-balance gravity. Gravitational collapse is one of the most 
important problems in general relativity. According to the singularity theo- 
rems [T] there exist spacetime singularities in generic gravitational collapse. 
It has been an interesting problem to determine the nature of spacetime sin- 
gularity. The cosmic censorship hypothesis (CCH) [2] says that singularities 
appearing in gravitational collapse are always clothed by the event horizon. 

The final fate of gravitational collapse of the massive star depends upon 
the choice of initial data and equation of state. Many efforts have been made 
to check its credibility but no final conclusion is drawn. For this purpose, 
Virbhadra et al. [3] introduced a new theoretical tool using the gravitational 
lensing phenomena. Also, Virbhadra and Ellis [I] studied the Schwarzschild 
black hole lensing and found that the relativistic images would confirm the 
Schwarzschild geometry close to the event horizon. The same authors jS] 
analyzed the gravitational lensing by a naked singularity and classified it as 
weak naked singularity and strong naked singularity. In a recent paper [6], 
Virbhadra used the gravitational lensing phenomena to find the improved 
form of the CCH. 

Oppenheimer and Snyder [7] studied dust collapse for the first time and 
showed that singularity is neither locally or globally naked. This means that 
they found black hole as a final fate of the dust collapse. Eardely and Smarr 
[8] found that inhomogeneous model undergoes to gravitational collapse by 
forming a singularity that can be either locally or globally naked. 

There has been a growing interest to study gravitational collapse in the 
presence of perfect fluid and other general physical form of the fluid. Misner 
and Sharp [9] extended the pioneer work for the perfect fluid. Vaidya [10] 
and Santos [11] used the idea of outgoing radiation of the collapsing body 
and also included the dissipation in the source by allowing the radial heat 
flow. Markovic and Shapiro [T2] generalized the pioneer work with positive 
cosmological constant. Lake [T3] extended it for both positive and nega- 
tive cosmological constant. Sharif and Ahmad [H]-[I7J extended spherically 
symmetric gravitational collapse with positive cosmological constant for per- 
fect fluid. The same authors [18] have also investigated plane symmetric 
gravitational collapse using junction conditions which has been extended to 
spherically symmetric gravitational collapse |19j . 
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The behavior of electromagnetic field in strong gravitational field has been 
the subject of interest for the researchers over the past decades. According 
to Thirukkanesh and Maharaj [20], the inclusion of electromagnetic field in 
gravitational collapse predicts that the gravitational attraction is counter- 
balanced by the Coulomb repulsive force along with the pressure gradient. 
Sharma et al. [2TJ have concluded that electromagnetic field affects the value 
of red-shift, luminosity and mass of the relativistic compact objects. Nath 
et al. [22] have studied the gravitational collapse of non-viscous, heat con- 
ducting fluid in the presence of electromagnetic field. They concluded that 
electromagnetic field reduces pressure and favors the formation of naked sin- 
gularity. 

Recently, we have studied the effect of electromagnetic field on the gravi- 
tational collapse by taking spherically symmetric spacetime as interior region 
and Reissner-Nordstrom as exterior region of the star [23]. The present article 
investigates the previous work by taking the Friedmann universe models in 
the interior of star. In order to preserve the generic properties of the Fried- 
mann universe models in the presence of electromagnetic field, we follow 
[231123 an d assume that electromagnetic field is weak relative to matter, i.e., 
if E 2 is the electromagnetic field contribution in the system then E 2 << p. 
The main objectives of this work are the following: 

• To study the physical interpretation of electromagnetic field and cos- 
mological constant on gravitational collapse in the Friedmann universe 
models. 

• To see the validity of CCH in this framework. 

The plan of the paper is as follows: In the next section, the junction 
conditions are given. We discuss the solution of the Einstein-Maxwell field 
equations in section 3. The apparent horizons and their physical significance 
are presented in section 4. Section 5 presents the singularity analysis. We 
conclude our discussion in the last section. 

The geometrized units (i.e., the gravitational constant G=l and speed of 
light in vacuum c = 1 so that M = and k = = 8n) are used. All the 
Latin and Greek indices vary from to 3, otherwise, it will be mentioned. 
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2 Junction Conditions 



We derive conditions for the smooth matching of two regions (interior and 
exterior of a star) on the surface of discontinuity. For this purpose, we assume 
that S be a timelike 3D hypersurface which divides two 4D manifolds V~ 
and V + respectively. The interior manifold is taken as the Friedmann model 

ds 2 _ = dt 2 - a(t) 2 [d X 2 - f(d6 2 + sm6 2 d(f)% (2.1) 

where fk(x) is defined as 



fU) = <X, k = 0, (2.2) 




k — 1,0, —1 correspond to closed, flat and open models respectively. x(0 < 
X < oo for open and closed but < x < 71 f° r na t) is t ne hyper-spherical 
angle and a(t) is the scale factor. Further, x is related to radial coordinate 
r as follows: r = sinx (closed), r = x (flat) and r = sinhx (open). The 
Reissner-Nordstrom de-Sitter spacetime is taken as the exterior manifold 

ds 2 + = ZdT 2 - -dR 2 - R 2 {dd 2 + sin 9 2 d(f) 2 ), (2.3) 

where 

M and A are constants and Q is the charge. 

The junction conditions are given as follows [26J: 

1. The continuity of first fundamental form over E gives 

(ds 2 _ ) s = (ds 2 + )x = ds 2 ,. (2-5) 

2. The continuity of second fundamental form (extrinsic curvature) over 
£ yields 

[Kij] = K+ - = 0, (ij = 0,2,3) (2.6) 
where Kij is the extrinsic curvature defined as 
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, +/ d 2 xl dx*idx v +s , , N 

Here £° = t, £ 2 = 9, £ 3 = are the corresponding parameters on E, x± stand 
for coordinates in V ± , the Christoffel symbols Y a are calculated from the 
interior or exterior spacetimes and are the components of outward unit 
normals to E in the coordinates x a ±. 

The equation of hypersurface in terms of interior spacetime V~ coordi- 
nates is 

/-(*,*) = X-Xe = 0, (2.8) 

where xs is a constant as E is a comoving surface forming the boundary 
of interior matter. Also, the equation of hypersurface in terms of exterior 
spacetime V + coordinates is given by 

f + (R,T) = R-R s (T) = 0. (2.9) 

When we make use of Eq. (12.81) in Eq. fl2.ll) . the metric on E takes the form 

(dsl ) s = dt 2 - a{tff(x^){d6 2 + sin# 2 d0 2 ). (2.10) 

Also, Eqs.(j23]) and Q yield 

1 d R 

(ds 2 + h = [Z(R S ) - ^y(^) 2 ]^ 2 - RUd6 2 + sinfl 2 # 2 ), (2.11) 
where we assume that 

z ^-zk0 )2> ° (2 - 12) 

so that T is a timelike coordinate. From Eqs.f l2.5p . (12. 101) and ( 12. lip , it 
follows that 

ife = (a/) E , (2.13) 

Also, from Eqs. (l2.8p and (I2.9p . the outward unit normals in V~ and V + , 
respectively, are given by 

n- = (0,a(t),0,0), (2.15) 

n+ = (-Rx,t, 0,0). (2.16) 



The components of extrinsic curvature Kf- become 



K 00 = 0, (2.17) 
K 22 = esc 2 OK^ = iff ah, (2.18) 

= (Af-TA-fgtf- + ±^ B> (2,9) 
iT+ = csc 2 6K+ = (ZRf)x, (2.20) 

where dot and prime mean differentiation with respect to t and x respectively. 
From Eq. fl2.6p . the continuity of extrinsic curvature gives 

K+ = 0, (2.21) 
K+ = K 22 . (2.22) 

Using Eqs. fl2Tn) - fT2~22l along with Eqs.(j22D, (12331) and (12141) . the junction 
conditions become 

(f% = 0, (2.23) 

Equations (1X13]) . fl2TT4j) . fT2T2l?j) and (EMD provide the necessary and suffi- 
cient conditions for the smooth matching of the two regions over E. 



3 Solution of the Einstein Field Equations 

In this section, we solve the Einstein field equations with cosmological con- 
stant for the Friedmann models containing the charged perfect fluid as the 
source of gravitation. The Einstein field equations with cosmological con- 
stant are given by 

G, u -Xg, u = K{T^ + T^). (3.1) 
The energy-momentum tensor for perfect fluid is 

Tp, = (p + p)Uf,u u - pg^, (3.2) 

where p is the energy density, p is the pressure and = 5® is the four- vector 

velocity in co-moving coordinates. TJ]^. is the energy-momentum tensor for 
the electromagnetic field given by 

T inT ] = -^{-9 5U1 F, 5 F VU3 + \g^F Suj F s "). (3.3) 
6 



With the help of Eqs.fl32J and flS^l . Eq.f l3~Tj) takes the form 



R» v = 8vr[(p + v)u,u v + i(p - p)^ + T£ m) - i^T^)] - A^. (3.4) 
Now we solve the Maxwell's field equations 

Fpv = 4> v ,n - (f>w, (3.5) 
F^v = AttJ^ (3.6) 

where M is the four potential and J M is the four current. Since the charge 
is at rest in this system, the magnetic field will be zero. Thus we can choose 
the four potential and four current as follows 

<^ = (0(t,r), 0,0,0), J» = au», (3.7) 

where a is charge density. Using Eqs.f l3.5p and ( 13. 7ft . the non-zero compo- 
nents of the field tensor are given as follows: 

Fqi = —F w = —-J-- (3.8) 
dx 

Also, from Eqs.fjMD, (D an d (133]) , we have 

+ 2y = W, (3.9) 

a ^ + % = a (310) 

Integration of Eq. fl3.lip implies that 

90 1 



dx a/ 1 



q(x), (3-11) 



where q(x) = An f* aa 3 f 2 dx, is the total charge distribution in the interior 
spacetime. This amount of charge is the consequence of law of conservation 
of charge, i.e., J!t = 0. It is clear that Eq. fl3.10p is identically satisfied by 
Eq. 113. lip . The electromagnetic field intensity is given by 

E = AV2- ( 3 - 12 ) 
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Equations f l3~TTj) and f[3~T2j) yield 

^ = aE. (3.13) 

ox 

Using Eqs. (El]) and fl3~T3l) . we get 

Foi = -i^io = -a^- (3.14) 
The non-zero components of T^J^ and its trace free form turn out to be 

rp{em) _ 1- 77>2 T (™) _ '- t H _ ^ tti2/_j-\2 

i 00 - > i ll - "^-^ fl ' ^22 - l a /J > 

<57T <57T 57T 

T 3 ( 3 em) = T 2 ( 2 em) sin 2 T (e-) =0 . 

When we use these values, the field equations (13. 4p for the interior spacetime 
takes the form 

R 00 = -3- = 4tt( P + 3p) + E 2 - A, (3.15) 

Ru = -- - + —J-j — ^{P - p) + E 2 - A, (3.16) 

R22 = - a --(y + ^ + £) 2 -h=Mp-p)-E 2 -A(3.U) 
a a a j j j 

R 33 = sm 2 6R 22 , (3.18) 

We would like to mention here that all the results are valid for E 2 « p 
and hence for stiff matter (p = p), E 2 << p. Integrating Eq. (12.231) with 
respect to t, it follows that 

/' = W, (3.19) 

where W = W(x) is an arbitrary function of \. The energy conservation 
equation 

t;. v = (3.20) 

for the perfect fluid with the interior metric shows that pressure is a function 
of t only, i.e., 

p = p{t). (3.21) 

Using the values of /' and p from Eqs. fl3Tl9l and fEOTD in Eqs. fl3TT5D -f l3T7D . 
it follows that 

a ,a, 9 (1 — W 2 ) o , . , 

2- + - 2 + 1 . — A + E — 8np(t). 3.22 

a a \ a J) 
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We consider p as a polynomial in t as given by 

p(t) = p c (|)~ s , (3.23) 

where T is the constant time introduced in the problem due to physical reason 
by re-scaling of t, p c and s are positive constants. Further, for simplicity, we 
take s = so that 

pit) = p c . (3.24) 

Now Eq. fl3T2"2"|) gives 

a (1 — W 2 ) , . 

2- + R 2 + , ^2 = A + E ~ 8lT Pf 3 " 25 

a a \ a J) 



For the static charges E is taken as time independent [27j . so integration of 
above equation with respect to t, yields 

FT 2 — 1 9 s a 2 m , 

a 2 = —pt- + (A + E 2 - 8np c )- + 2—, (3.26) 
J 3 aj 

where m = m(x) is an arbitrary function of x an d is related to the mass of 
the collapsing system. Substituting Eqs.f l3.llJ]) . (I3.26P into Eq. (13. 15j) . we get 

2E'E 

m' = ^-(«/) 3 + a 3 ff 2 [Air(p c + p) + 2E 2 ]. (3.27) 

For physical reasons, we assume that (p c + p) ^ 0. Integrating Eq. (13.271) 
with respect to x, we obtain 

m( X )=47ra 3 / (p + Pc )f'f 2 d X + 2 / E 2 f'f 2 d X +-a 3 E'Ef 3 d X + m , 
Jo Jo d Jo 

(3.28) 

where mo is taken equal to zero because of finite distribution of matter at 
the origin. The function m(x) must be positive because m(x) < implies 
negative mass which is not physical. Using Eqs.f l3.19l and (13. 26ft into the 
junction condition Eq. (l2.24p . it follows that 

M=-^j + m+^(A + E 2 - 8np c )(af) 3 . (3.29) 



9 



The total energy M(x, t) at time t inside the hypersurface S can be evaluated 
by using the definition of mass function with the contribution of electromag- 
netic field for the Friedmann model, which is given by 

M( X ,t) = \{af){l + (df) 2 - f) + £j. (3.30) 

Replacing Eqs. fl3TT9|) and f l3T26|) in Eq. fl3T30|) . we obtain 

M(r, t) = m{r) + (A + E 2 - 8rcp c )^- + (3.31) 

6 2a/ 

From Eqs. fl3T29|) and fETBTj) . it can be found that M(r,t) = E M if and only 
if q = Q. This result provides the necessary and sufficient conditions for the 
continuity of mass in the interior and exterior regions over boundary surface 
E. 

Now we take (A + E 2 — 8np c ) > such that E 2 « 8-np and assume that 

W( X ) = 1. (3.32) 

In order to obtain the analytic solutions in closed form, we use Eqs. fl3.19p . 
fl3T26|) and fl332|) so that 

(af) = ( A ° o )^inha«( Xi t) (3.33) 
A + E z — 8np c 

where 



. i •, V X 3(A + E 2 — 8irp c ) 
a{x,t) = [t s (x)-t)\. (3.34) 

Here t s (x) is an arbitrary function of x an d is related to the time of formation 
of singularity. 



4 Apparent Horizons 

In this section, we discuss the formation of apparent horizons. The boundary 
of two trapped spheres whose outward normals are null is used to find the 
apparent horizons. Moreover, we discuss the the physical significance of 
apparent horizons 1.6., cLFGcL of apparent horizons, time difference between 
apparent horizons and singularity etc. For the interior spacetime, we find 
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the boundary of two trapped spheres whose outward normals are null as 
follows: 

<r{*f)A*f)*> = (of) 2 - U'f = o. (4.1) 

Using Eqs. fl3.19p and (13.261) in this equation, we get 

(A + E 2 - %ixp c )(aff - 3(o/) + 6m = 0. (4.2) 
When A = 8irp c — E 2 , we have (af) = 2m. This is called Schwarzschild 

3 



horizon. For m = p c = K = 0, we have (af) = y -|, which is called de-Sitter 
horizon. Equation ( 14. 2 p can have the following positive roots. 

Case (i): For 3m < , 1 =. we obtain two horizons 

(a/)c = 7Wrw=m m l (4 - 3) 

where 



cosy? = -3m v /(A + E 2 -8vrp c ). (4.5) 



If we take m — 0, it follows from Eqs. (l4.3p and (I4.4p that (af) c — y ( A+E 2_ 8wp ) 

and (af)b = 0. (af) c and (af)b are called cosmological horizon and black 
hole horizons respectively. For m^0 and A ^ 8Trp c — E 2 , (af) c and (af)t, 
can be generalized (28] respectively. 

Case (ii): For 3m = , 1 , there is only one positive root which 
corresponds to a single horizon i.e., 

(o/)c = (af) b = t= = = = (a/)*. (4.6) 
V(A + £ 2 - 8vrp c ) 

This shows that both horizons coincide. The range for the cosmological and 
black hole horizon can be written as follows 



< (af) b < 1 = < (af) c < J ; 7 . ( 1.7) 

a/(A + £ 2 - 8vrp c ) V (A + £ 2 - 87rp c ' 
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The black hole horizon has its largest proper area 4n(af) 2 = ( A+E 2^ 87rp \ and 
cosmo logical horizon has its area between ,. . J? a — T and ,. . }% n a — r. 

Case (iii): For 3m > , 1 , there are no positive roots and con- 

V ' v /( A+£ ;2_ 87rpc) ' V 

sequently there are no apparent horizons. 

We now calculate the time of formation of the apparent horizon using 
Eqs.fllLM} and fl42|) 

t n = t s -—= 2 -sinh-V^^-l)^ (n = l,2). (4.8) 

v/3(A + E 2 — 8n Pc ) 1 2m ; 1 ; K J 



This implies that 

(fl/)n -cosh 2 a„, (4.9) 



2m 



where a n (r,x) = \/ 3 ( A+ ^ — [^s(x) — tn)]- Equations (I3.33H and (I4.6f) give 
( a f)c > ( a /)fc an d tb > respectively. The inequality t b > t c indicates 
that the cosmological horizon forms earlier than the black hole horizon. This 
condition confirms the formation of black hole. 

The time difference between the formation of cosmological horizon and 
singularity and the formation of black hole horizon and singularity can be 
found as follows. Using Eqs. (14 .3 [) - (14. 5 j) . it follows that 

d(^) 1 , sin | 3cosf s , 

- 3 + ? - \ < 0, (4.10) 



dm m siny cosy 

d (Mk) i sin (^) 3cos^^ 



dm m siny cosy 



>0. (4.11) 



The time difference between the formation of singularity and apparent hori- 
zons is 

r n = t s -t n . (4.12) 

It follows from Eq. (J3~9~j) that 



n 



^(lm) sinh a n cosh a n a/3 ( A + E 2 — &np c 



(4.13) 
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Using Eqs.f H~TU|) and (147131 . we get 



2m 



dm d(^^) dm m a/3(A + E 2 - 8np c ) sinh a% cosh a% 

. sin f 3 cos % , . 

x ^ + < 0. 4.14 

sm ip cos (/9 

This means that time interval between the formation of cosmological hori- 
zon and singularity is decreased with the increase of mass. Similarly, from 
Eqs.flHTJ andfllTSD, we get 

dr 2 1 



dm m a/3 (A + E 2 — 8rcp c ) sinh a 2 cosh a 2 

sin (^ 3cos^ x , x 

x ^— + ^— > 0. 4.15 

sm If COS </? 

This indicates that time difference between the formation of black hole hori- 
zon and singularity is increased with the increase of mass. 



5 Singularity Analysis 

The Riemann tensor is used to determine whether a singularity is essential 
or removable. If the curvature becomes infinite at certain point, then the 
singularity will be essential otherwise removable. Many scalars can be con- 
structed from the Riemann tensor but symmetry assumption can be used to 
find only a finite number of independent scalars. Some of these are 

p p „ab p p p Dab p p ryabcd p pat pcd 

ti\ — ti — g ti a b, ti 2 — tiabti , -Tt3 — tiabcdti , ti 4 — K cd ti ab . 

Here, we give the analysis for the first invariant commonly known as the 
Ricci scalar. For the Friedmann model, it is given as 

p ~3aa/ 2 + 2f"f - 3a 2 / - 1 + f 2 

R= — . 5.16 

a 2 f 

By definition a > and - > [2H], it follows that curves of a(t) versus t 
must be concave downward and must reach a(t) = at some finite time in 
the past. Let us recall this time t = at which R = oo. In cosmology, 
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extrapolation of the universe expansion backwards in time yields an infinite 
density at finite past. Also, if the strong energy condition [TJ is satisfied, i.e., 
p + p > and (p + 3p) > then a = at t = which implies the divergence 
of scalar curvature polynomial where p — > oo. This is spacelike singularity 
usually called big bang singularity or initial singularity [30] . 

6 Summary and Conclusion 

In this paper, we have analyzed the charged perfect fluid collapse with pos- 
itive cosmological constant in the Friedmann models. For this purpose, 
we have found junction conditions between the Friedmann models and the 
Reissner- Nordstrom de-Sitter spacetime. The junction conditions provide the 
gate way for the exact solution of the field equations with interior spacetime 
(Friedmann models). The solution of the field equations helps to discuss the 
dynamics of the collapsing system as follows: 

The acceleration parameter a /a, given by Eq. fl2.14j) . will be zero, positive 
or negative for 4n(p + 3p) + E 2 - A = 0, 47r(p + 3p) + E 2 - A < or 
47r(p + 3p) + E 2 — A > respectively. The variation of the scale factor a{t) 
helps to describe the different stages of matter in the Friedmann models of 
the universe [31]. If the scale factor a{t) is decreasing, i.e., a{t) < there will 
be collapsing (contracting) phase. For increasing scale factor i.e., a{t) > 
we have the expanding phase while the point where a(t) =0 corresponds to 
bounce point. Consequently, the Hubble parameter will be H < 0, H > 
and H = for collapsing, expanding and bouncing phases respectively. Also, 
we can conclude the following: 

• The Newtonian force and acceleration of matter have the same value 
over the hypersurface S, i.e., (~™ + (A + E 2 - 87rp c )^) s (see [23] 
for detail). In this case, the repulsive force can only be generated 
if A > (87rp c — E 2 ) such that 87rp c » E 2 over the entire range of 
the collapsing sphere. In the case of charged perfect fluid collapse 
with Tolman-Bondi spacetime [23] there is no restriction on matter and 
electromagnetic field then the results are valid only for A > (87rp c — E 2 ) 
such that 87rp c > E 2 . It is clear that in the first case the cosmological 
constant attains higher value than the later case. Thus the cosmological 
constant plays an effective role to slow down the collapse in the present 
case than previous one. In other words, isotropy and homogeneity of 
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matter causes to introduce resistance against collapse in the presence 
of charge. 

• Since the cosmological constant A is affected by pressure and electro- 
magnetic field, we can say that electromagnetic field reduces the effects 
of A as compared to perfect fluid case by putting the restriction on A. 
Hence electromagnetic field increases the gravitational collapse as it 
decreases the repulsive force produced by A. 

• Two physical horizons (cosmological and black hole horizons) are found 
whose area is decreased by cosmological constant and electromagnetic 
field. It follows from Eq. fl4.8p that both horizons form earlier than 
singularity, so singularity is covered (back hole) and CCH seems to be 
valid in this case. 

• Time difference between the formation of apparent horizon and sin- 
gularity is decreased by electromagnetic field. Thus we can say that 
singularity must form earlier than the apparent horizons. Hence elec- 
tromagnetic field favors the formation of naked singularity. But such 
situation can never occur because electromagnetic field does not play 
the dominant role in this case. 

• It is found that the time difference between the formation of cosmolog- 
ical (black) horizon and singularity is decreasing (increasing) function 
of mass of the collapsing system. 
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